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ABSTRACT 

We identify a quantity in the c = 1 matrix model which describes the wavefunc- 
tion for physical scattering of a tachyon from a black hole of the two dimensional 
critical string theory. At the semiclassical level this quantity corresponds to the 
usual picture of a wave coming in from infinity, part of which enters the black hole 
becoming singular at the singularity, while the rest is scattered back to infinity, 
with nothing emerging from the whitehole. We find, however, that the exact non- 
perturbative wavefunction is nonsingular at the singularity and appears to end up 
in the asymptotic region "behind" the singularity. 
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In a previous paper it was suggested that the one dimensional matrix model 
may be used to study the nonperturbative behavior of massless "tachyons" cou- 
pled to the black hole background of the two dimensional critical string' 3 " 41 . This 
follows from the fact that at the semiclassical level a certain integral transform of 
the fluctuation of the collective field of the matrix model M 161 around the classical 
ground state value satisfies the same linearized equation of motion as that of the 
massless tachyon in a blackhole background. Indeed this connection between the 
matrix model and the black hole is essentially the same as that between the liou- 
ville theory and the SL(2, R)/U(l) coset model discovered earlier by Martinec and 

ryi rgi 

Shatashvilli . It was soon found in that a similar transformation exists in the 
bosonized form of the fermionic field theory of the matrix model developed in 
and leads to the same black hole interpretation. A large class of transformations 
between collective field theory and black hole wavefunctions have been discussed 
in 1101 . In fact, there are several pieces of evidence from continuum treatments that 
the liouville and black hole backgrounds are closely related 1111 . 

The presence of this connection does not imply that the physics of the black 
hole background is identical to that of the liouville background. Rather, at the 
semiclassical level, the integral transform allows one to obtain the behaviour of a 
tachyon in a black hole background from matrix model quantities just as singular 
gauge transformations in gauge theories allow one to obtain wavefunctions in pres- 
ence of vortices or monopoles from those in free space. However, if one examines 
the structure of the integral transform, one finds that near the singularity of the 
black hole this receives contributions from the collective field in regions where the 
coupling is strong. Far from the black hole the semiclassical approximation is good 
so that one can meaningfully talk about incoming and outgoing tachyon states. 

Beyond the semiclassical expansion the matrix model quantities transformed as 
above would represent correlations in some quantum field theory of scalar particles 
whose semiclassical behaviour is identical to that of tachyons moving in a static 
black hole. We could thus use matrix model results to study non-perturbative 
behaviour of tachyon propagation. Indeed one feature of such non-perturbative 
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effects have been observed in where it was shown that, with the particular defi- 
nition of this transform used in that paper, the classical value of the background is 
singular at the location of the black hole singularity, but the exact value is com- 
pletely regular. A similar behavior has been found for the transform of a class of 

[12] 

wavepackets which represent "spikes" on the fermi sea . However these spikes 
transform into tachyon wavefunctions which do not represent scattering processes 
as described above. 

What is lacking in all the papers quoted above is a clear identification of 
the quantities which have to be computed in the matrix model which represent 
a physical scattering process in the black hole background. For example, at the 
semiclassical level, the transforms defined in either [1] or 181 represent a process where 
there is a non-zero flux emerging from the past null infinity and hence do not 
represent a process in which the scalar particle comes in from infinity, part of the 
wave being scattered and part of the wave absorbed by the black hole. One of the 
transforms considered in' 101 does represent a scattering type wavefunction in the 
exterior region, but not in the interior region - where the wave has a component 
which does not vanish on the extension of the past horizon. 

In this paper we make this identification by noting that the large class of 
transforms which take matrix model quantities into single particle wavefunctions 
in a black hole background correspond to different boundary conditions' 7 " 101 . We 
find the correct linear combination of these transforms which corresponds to a wave 
vanishing on the past horizon and hence represents a physical scattering process 
from a black hole. We show that, as expected, single particle wave functions which 
represent such a scattering process diverge at the black hole singularity at the 
semiclassical level. We then consider quantum corrections to this single particle 
wave function using the exact results of the matrix model. We show that while to 
all orders of string perturbation theory the divergence of the wavefunction at the 
location of the singularity persists, the exact non-perturbative answer is completely 
nonsingular. We show that in the aymptotic region "beyond" the singularity this 
wave corresponds to a purely outgoing wave. Thus at the exact quantum level the 
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particle thrown in from the standard asymptotic region in the exterior of the black 
hole ends up in a different asymptotic region beyond the black hole. 

The semiclassical limit 

Consider the two types of loop operators in the matrix model of a matrix M y - (t) 
The macroscopic loop operator with real loop lengths is defined by the quantity 

oo 

W(p,t) = Tr e~ pM = J dx e~ px (j>{x, t) (1) 

We will also define macroscopic loop operators with imaginary loop lengths which 
are defined by 

oo 

V(z, t) ee Tr e izM = J dx e izx f) (2) 

— oo 

The field (/>(x,t) is the density of eigenvalues of M y -, the collective field. In 

[13] 

terms of the fermionic fields ip(x,t) of the fermionic field theory of one has 
4>(x,t) = ip*(x, t)ip(x, t). In the fermionic field theory the correlations of W(p,t) 
may be obtained by first computing the correlations of V(z, t) and then performing 
a suitable analytic continuation, as explained in [14] . However, in the considerations 
which follow it is important to keep in mind the difference between the two quan- 
tities. 

In the lowest order of the semiclassical expansion one can obtain linearized 
equations for the matrix elements of the loop operators between one particle states 
and the ground state, which we will denote by V and W. Semiclassically these 
may be obtained by expanding the collective field <f>(x, t) around its classical ground 
state value 4>o(x,t) = ^(x 2 — 4/x)a 

(f)(x,t) = (f) (x,t) + d x r](x,t) (3) 

and replace <f>(x,t) in (1) and (2) by d x rj(x,t). Using the equations of collective 
field theory it may be seen that at the classical level the linearized equation of 
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motion for r)(x, t) leads to the following linearized equation for W(p, t) and V(z, t) 



[( P d p y-d{-4^}w(p,t) = o 

( 4 ) 

[(zd z ) 2 -d 2 + 4fiz 2 }V(z,t) = 



The first equation in (4) is identical to the Wheeler-de Witt equation obtained 
in the liouville theory. Indeed using the semiclassical expression for energy eigen- 
states r)(x,t) = e lu}t sin(co cosh^ 1 ^-^/^)) °f the matrix model satisfying the correct 
Dirichlet boundary condition the fluctuation of W{p,t) is given by the modified 
Bessel function e Ki v {2^JJl p) which represents (in the coordinate log p) a wave 
coming in from infinity and getting reflected perfectly from the "liouville barrier" . 
This is of course the basic picture of scattering in the liouville theory. Formally 
V(z,t) may be regarded as the liouville wave function with negative cosmological 
constant. We shall not, however, assign any physical meaning to V(z,t) and we 
will not need any. 

Consider now the following classes of transforms of V and W 

oo oo 



T±(u,v) = J dz J dt e ±>*(e'«-e~'«) V(±z,t) 

-oo 
00 oo 

S±(u,v) = J dp J dt e ±^v+e-*u) w{p,t) 



(5) 



-oo 



Using equations (4) it may be easily verified that each of these quantities T±, S± 
satisfy the following linearized equation 

[4(uu + fx)d u d v + 2(ud u + vd v ) + l]Xi(u, v) = (6) 

where X\ stands for any of the T±, S± This is precisely the linearized equation of 
motion for the massless tachyon of the two dimensional critical string moving in a 
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black hole background. The coordinates (u, v) are the Kruskal-like coordinates in 
terms of which the background metric and dilaton fields are* 



G U v — G vu — 2 x G uu — G vv — 

uv + a) (?) 

12 

D(u, v) = -- log ( — uv + a) 
2 a 

The black hole mass is given by a = ^u. In terms of the coordinates (w, v) the 
various regions of the black hole geometry are 



u — re v — r e for u,v >0 Region I 

ft Or 

u = —r e v = r e for w < 0, t> > Region II 

u = —r e v = — r e for u, v < Region III 

u = r e 9 v = — r e~ 9 for u > 0, v < Region IV 



(8) 



Thus the Region II contains the future singularity while the Region IV contains 
the past singularity. 

The transform defined in' 11 is given by S+, while the transform defined in 181 is 
a specific linear combination of all of these independent transforms^. Given a 
space-time history of <f>(x, t) these various transforms provide different space time 
histories of a tachyon in a black hole background. Semi classically these histories 
may be easily evaluated using the expression for the fluctuations of the collective 
field given above and are given by hypergeometric functions, as done for S + m [1] . 
However, the different transforms lead to different hypergeometric functions, i.e. 
waves with different boundary conditions' 101 . For example S+ with W chosen to 

* Our conventions are those of ' 

f In 181 the transform involves a collective field which is the density of fermions in the momen- 
tum space conjugate to the eigenvalue coordinate. However, in the double scaling limit the 
interchange of momenta and coordinates simply changes the sign of the cosmological con- 

rgj 

stant. This may be then used to relate the part of the transform of involving momentum 
space collective fields in terms of the quantities we have defined above. 



6 



be of the form W ~ e is given in Region I by 



S + (r, 6) =—&A(v)F(- - iu, -; 1 - «/; --)+ 

if Zi Zi IX 



(9) 



where A(z/) = and we have omitted an overall //-dependent constant. It 

is clear from (9) that near the horizon at r = this represents a linear combination 
of a wave coming out of the past event horizon as well as one going into the future 
event horizon. This is not what we want to study for the scattering of a particle 
from a black hole. For the latter physical problem we must ensure that there 
is nothing emerging from the past horizon. We specifically want to study what 

[8] 

happens to the wave at the singularity. The transform of also consists of a wave 
which is nonvanishing at the past horizon. 

It is indeed possible to write down a transform which represents a physical 
scattering of a tachyon from the black hole, i.e. a wave coming in from infinity - 
part of which gets inside the horizon while a part is scattered back with nothing 
emerging from the past horizon. This is obtained by first using the time translation 
invariance of the theory to work in terms of fourier transforms of the operators 
V(z, t) defined by 

V(z,t) = [ dv e ivt V{z,v) (10) 



The correct transform is then given, in terms of the fourier transformed quantities 
V(z,v) by 



oo oo 



T^U.V) = J dz J dt [e^ e ^(e t v-e- t u) + e - ¥ e -iz(e*«-e-*«)] jut ^ 
-oo 

(11) 

The above definition may be regarded as an operator definition. To extract single 
particle wave functions and hence obtain T sca t one has to take suitable matrix 
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elements which amounts to the replacement of V(z,t) by V(z,t). The integration 
over t in (11) may be explicitly performed separately in the different regions, with 
the result 

oo 

T^lt = e™ 6 J dz [K iu (2zr)][V(z, v) + V(-z, v)] Region I 
o 

oo 

f^ t = e M J dz [H\l\2zr) - H%\2zr)][V(z, v) + V(-z, v)\ Region II 
o 

(12) 

The expressions in the other regions may be trivially obtained from these by chang- 
ing the sign of one of the coordinates appropriately. 

Let us first evaluate this transform semiclassically using collective field theory 
as discussed above. This leads to 

V sc (z,v) = A(v) H\l\2^Jlz) V sc (-z,v) = -A(v) H^{2^z) (13) 

where H\J denote the standard Hankel functions, and A{y) is a constant which 
is irrelevant to our discussion. The superscript sc stands for semiclassical. Substi- 
tuting these expressions in (12) and performing the integral over z lW] one obtains 

v) = B(u) v ~™ F{\ - iv, I; 1 - iv- -^) (14) 

where B(v) = — t ^ r ^^^_ lS ^ /l ^ 7r ' y ' ) , and F stands for the hypergeometric function. 
This expression is valid throughout Region I while in Region II it is valid for 
— uv < fi (note in Region II uv is negative), i.e. in the region between the horizon 
and the singularity. In the exterior region this represents a wave which vanishes 
on the past horizon and may be written (by using standard relations between 
hypergeometric functions) in terms of a specific linear combination of right and 
left moving waves at an infinite distance away from the horizon. Thus we have an 
incoming wave and an outgoing wave, but nothing emerging from the past horizon 
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- precisely what we wanted. At the black hole singularity uv = —/j, the wave is 
singular. The singularity comes from the upper limit of integration in (12). For 
large z the asymptotic forms of the Hankel functions are 



H$ ~ ±= ^ ¥U) Ht ] ~ 4= e~^~^ (15) 

W yfz \/Z 

Using this it may be seen that in (12) the integral of products H^(2rz) 2^/Jiz) 

(2) (2) 

and H\ v ' {2zr-)H\J {2^fjiz) are finite while the cross terms lead to a logarithmic sin- 
gularity at r = yjji. Therefore in the semiclassical limit a part of the wave that 
came in from infinity and went inside the horizon "crashes" into the singularity*. 

The Exact wavef unctions 



For non-perturbative computations of V(z, t) we use the fermionic field theory 
of the matrix model' 131 1141 . As explained in 1141 , this is done as follows. Consider the 
two point function < V(z\, v)V(z2, —v) >. The single particle tachyon wavefunc- 
tion in the standard interpretation of the matrix model is obtaining by computing 
the quantity 

V(z 2 , ")=<f^ H iV ~ X < V(z!, u)V(z 2 , -v) > (16) 

and then performing a suitable analytic continuation to real loop lengths. The 
reasoning behind this identification is that the various powers of z\ or Z2 denote 
contributions from various operators and the above contour intgeral simply isolates 
the part which is the contribution of the tachyon operator as dictated by the 
correct scaling dimension. (We will restrict our attention to nonintegral v. The 
interpretation for integer v is more involved, see e _g_ [17I18I19] anc [ re i a ted to the 
symmetry structure of the theory' 201 [21] 1181 [22) ) It may be also simply verified that 

* The transform considered in 111 does not diverge at the singularity. However, as noted above 
this does not describe a physical scattering process with the correct boundary conditions. 
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the quantity V(z, v) satisfies the Wheeler de Witt equation with negative /x at the 
semiclassical level. Thus after analytic continuation z — > il (in the sense explained 
in Ref[14]) this does represent the physical tachyon single particle wave function. 

We are, however, not interested in extracting liouville wave functions; we will, 

+ 

therefore, make no such analytic continuation 

The two point function in (16), or rather its /i-derivative may be computed 
exactly in the fermionic field theory, either directly as in [141 or using Ward identities 
to relate it to the one point function as shown in' 191 . The result is (after a change 
of integration variables in the formulae of these papers) 

oo 

#- < V(z l7 u)V(z 2 , -v) >= Im / , §— r e *+i coth «A*>x 

o/x J fj,sinh{£/2fj,) 

o 

™ sinh(v€/2(x) z\z 2 , ^ Ai r r sinh{r£/2) z\z 2 „ 

1 sin(iru) iv[ sinh(Z/2ixy + r 2 + z/ 2 r[ sinh^/2ix) n 

(17) 

The semiclassical expansion is obtained by expanding in powers of the string cou- 
pling i keeping I = y/jlz fixed. From (17) it is clear that this means replacing the 
hyperbolic functions by their power series expansions. To the leading order it is 
easy to see that the answers depend on the combinations yfJLz\ and y/jxz2, as one 
would expect on the basis of the Wheeler de Witt equations. 

The single particle wave functions V(z,u) comes from the first term of the 
series expansion of the Bessel function Ji v (z\Z2 / sinh^ / 2)) , which goes as [z\z-i) w . 
The result is 

d,V(z,u) = Im A{y) {z/2) w j — ( ^ 2)]iy+1 coth |] smhH /2) (18) 



Try 

where A{y) = ^j^^^f^/j ■ The integration over £ may be performed in terms 



f We are using a definition of the nonperturbative theory which has a potential behaving as 
—x 2 both for positive and negative x 1231 . Since we do not make any analytic continuation, 
the quantities we define are unambiguous as opposed to non-perturbative S-matrices which 
have to be defined by analytic continuation to real loop lengths. 
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of Whittaker functions . The result is 

d,v(z, v) = im B(i/)i[r(i-vi) w^Qto (-^ 2 )-r(^-i//+i^)^ (M _ |)if (-^ 2 )]| 

(19) 

where Vl^ K stands for the Whittaker function. For zero energy v — this is the 
derivative of the vacuum expectation value of the operator V(z,t) as expected. 
This takes the simple form 

V(z, 0) = Im [e^ T(i - iA*) - WWH* 2 )] (20) 

Several properties of these exact expressions for V(z, v) are worth noting and 
will be useful in a moment. First, to take the semiclassical limit one has to remem- 
ber that I = ^[Jiz has to be kept fixed while ^//Z — > oo. Using standard asymptotic 
forms of Whittaker functions 1241 it may be checked that the above expression reduce 
to the semiclassical expressions involving Hankel functions in this limit. Secondly, 
for small values of z these wavefunctions reduce to the semiclassical wavef unctions 
in the corresponding limit. This may be seen from (18) from where it is clear that 
for small z only large values of coth |, or only small values of £ contribute. In 
that case one can replace the hyperbolic functions by their power series expan- 
sions. But this is precisely what yields the semiclassical limit. Alternatively one 
can directly take the result (19) and use the small- z behaviour of the Whittacker 
functions 1241 and verify that one indeed recovers the same answer as one obtains 
from the small z behaviour of the Hankel functions involved in the semiclassical 
wave function. This is simply a reflection of the fact that the coupling constant 
in the string field theory of the matrix model goes rapidly to zero away from the 
wall. In our case this region corresponds to the region of small z. 

The quantity d^T sca t(u, v) can be now calculated by inserting (19) in (12). The 
first point to note about the reaulting expressions for the exact T sca t(u,v) is that 
for large values of \uv\, i.e. for large values of r the exact expressions agree with 
the semiclassical expressions in this region. This follows from the large r behavior 
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of the Bessel functions in (12). However, we have just argued that for small z 
the exact wavefunctions V(z, v) approach their semiclassical estimates. Thus for 
large r the black hole wavefunctions also agree with their semiclassical answers. 
This fact shows that in the asymptotic region of the black hole the coupling of the 
theory is indeed weak. By the same token we expect that near the black hole there 
will be significant deviations from the semiclassical behavior. 

We saw in that in the semiclassical limit the logarithmic singualrity of the 
scattering wavefunction T sca t came from the upper limit of integration, i.e. from 
large z. This came about from a competition between the exponentials coming 
from the matrix model wavefunctions and the Bessel functions in (12). However it 
is precisely in the large z region that the matrix model wavefunctions differ most 
from the semiclassical answers. In fact, the asymptotic behavior of the Whittaker 
function is given by 

W X:K (-iz 2 ) ~ z 2X (21) 

and there is no competition between this and the e 2lzr coming from the Hankel 
functions. The result is in fact finite at the singularity. 

It is important to note that the scattering wavefunction remains infinite at the 
location of the black hole singularity to all orders in the semiclassical expansion. 

[81 

This is similar to the result of about the value of the one point function. The 
smearing of the behavior of the wavefunction near the singularity is a genuinely 
nonperturbative effect. 

Discussion 

What does this result mean ? Our matrix model black hole connection shows 
that one can obtain wavefunctions (and presumably the S-matrix) of tachyons 
moving in a black hole background at the semiclassical level by evaluating suit- 
able quantities in the matrix model. Given this connection we evaluate the same 
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quantities non-perturbatively in the matrix model. Our interpretation of these 
non-perturbative answers is that we are evaluating correlations or wavefunctions 
in some (possibly nonlocal) theory of a scalar field T(u, v) whose semiclassical 
behavior is identical to that of tachyons in a black hole background. 

We found that while the exact wavefunctions approach the semiclassical wave- 
functions in the asymptotic region, they depart significantly inside the black hole, 
which is therefore a region of strong coupling. In this region it is not very mean- 
ingful to expand around the standard classical ground state solution and interpret 
the resulting fluctuations as "particles". Nevertheless T sca t(u,v) can be correctly 
interpreted as single particle tachyon wavefunctions in the aymptotic regions and 
hence it is meaningful to talk of a scattering process. 

As emphasized in [lol the region "beyond the singularity", i.e. the portion of 
Region II with — uv > \x has to be taken seriously for the string theory black 
hole. Since the exact answer for T sca t(u,v) is nonsingular it is now meaningful to 
continue it across the singularity and ask what happens "beyond the singularity" . 
For large r in Region II we can either look at the small z behavior of the Whittaker 
functions' 241 or simply look at the semiclassical answer. This can be easily computed 
(the answer (14) is not valid here since — uv > jj) 



which shows that the wavefunction has support on one of the null infinities (corre- 
sponding to large values of v). The wave simply lands up in the asymptotic region 
beyond the singularity. 

To really see whether we are describing the behavior of tachyons of two dimen- 
sional string theory in a black hole background, one has to compute S'-matrices 
and compare them with continuum S'-matrices, e.g. those obtained from the coset 
model. After all, the identification of quantities in the matrix model with tachyon 
fields in a liouville background can be considered as correct because the tree level 
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S-matrix computed by matrix model methods ° agreed with those obtained from 
the continuum' 261 . Without a similar continuum calculation in a black hole back- 
ground the relationship between the theory of tachyon fields which result from the 
transform discussed above and string theory is unclear. Nevertheless our results 
show that in this exactly solvable theory of tachyons nonperturbative effects may 
drastically alter the behavior near the singularity. 

In the standard interpretation of the matrix model it is troublesome to assign 
a clear spacetime interpretation in the strongly coupled region 191 [14] . For example 
there is no relativistic invariance in this region. In our case this implies that the 
region inside the black hole does not have a clear relativistic space time interpre- 
tation. This is, in fact, expected in a string theory. As noted above, the fact that 
the theory is weakly coupled far from the black hole allows us to meaningfully 
desribe scattering processes involving tachyons. However this makes the absence 
of singularities of the exact wavefunction difficult to understand. Assuming that 
our complicated field theory of massless "tachyons" is indeed a string theory, the 
only statement we can make is that the string does not "see" any singularity even 
though it perceives a singularity at the semiclassical level. Since the graviton and 
dilaton fields are not explicitly present in the formalism it is confusing to ask 
whether the singularity itself has vanished. 

Finally, it would be very interesting if one can understand time dependent 
backgrounds in string theory, like formation of black holes, in terms of the matrix 
model. Since there is no explicit field for the graviton or the dilaton it would not 
be possible to see the black hole "evaporate away" , but its effects on the scattering 
of tachyons can still be studied. 
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